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Abstract 

In three dimensions, the construction of bi-Hamihonian structure can be 
reduced to the solutions of a Riccati equation with the arclength coordinate 
of a Frenet-Serret frame being the independent variable. Explicit integration 
of conserved quantities are connected with the coefficients of Riccati equation 
which are elements of the third cohomology class. All explicitly constructed 
examples of bi-Hamiltonian systems are exhausted when this class along with 
the first one vanishes. The latter condition provides integrating factor for ex- 
plicit integration of Hamiltonian functions. For the Darboux-Halphen system, 
the Godbillon-Vey invariant is shown to arise as obstruction to integrability of 
integrating factor. 



1 Introduction 

A Poisson structure on a manifold is defined by a skew symmetric contravariant 
bilinear form subjected to the Jacobi identity expressed as the vanishing of the 
Schouten bracket of Poisson tensor with itself This structure having no 

non-degeneracy requirement becomes the basic underlying geometry to study 
non-canonical Hamilton's equations on odd dimensional manifolds as well as the 
Hamiltonian structures of nonlinear evolution equations [3]- [5]. 

The first interesting case of a completely degenerate finite dimensional Hamil- 
tonian structure occurs in three dimensions. Many works have been devoted 
to the study of three dimensional dynamical systems with primary concern on 
quantization, construction of conserved quantities, Hamiltonian structures, inte- 
grability problems and their numerical integration using techniques from various 
areas such as Poisson geometry, differential equations, Frobenius integrability 
theorem and theory of foliations [6]-|26j. 



^ To Prof. Dr. Yavuz Nutku. The existence of local bi-Hamiltonian structures was somehow 
known to him while we were working on the Darboux-Halphen system. I remember asking 
him, to be sure not misunderstood what we have discussed, "you mean all these systems are 
bi-Hamiltonian?" He said, "of course they are". 
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In [37], we reduced the problem of constructing Hamiltonian structures in 
three dimensions to the solutions of a Riccati equation in moving coordinates of 
Frenet-Serret frame. All known examples of dynamical systems having two com- 
patible and explicit Hamiltonian structures are exhausted by constant solution. 
We concluded that in three dimensions vector fields which are not eigenvectors 
of the curl operator are at least locally bi-Hamiltonian. 

In this work, we shall extend the discussion to analysis of obstructions to the 
construction of global Hamiltonian structures. We shall present the Darboux- 
Halphen system as an example for which the Godbillon-Vey three-form can be 
explicitly worked out as an obstruction to the existence of global structure. 

2 Hamiltonian Systems in Three Dimen- 
sions 

We shall summarize the necessary ingredients of the bi-Hamiltonian formalism 
in three dimensions. See [B]-[3B] for details and examples. For x = {a;*} = 
{x,y,z) e R^, < g M and overdot denoting the derivative with respect to t, we 
consider the system of autonomous differential equations 

x-v(x) (1) 

associated with a three-dimensional smooth vector field v. Eq.([T]) is said to be 
Hamiltonian if the right hand side can be written as v (x) — (x) {dH (x)) 
where H (x) is the Hamiltonian function and ft (x) is the Poisson bi- vector 
(i.e. a skew-symmetric, contravariant two-tensor) subjected to the Jacobi iden- 
tity [fi (x) , f2 (x)] = defined by the Schouten bracket [1,. In coordinates, if 
di = d/dx^, the Poisson bi-vector is ft (x) = 57^'^ (x) dj A dk, with summation 
over repeated indices, and the Jacobi identity reads fl^^^dift'^^^ = where [jkl] 
denotes the antisymmetrization over three indices. It follows that in three di- 
mensions the Jacobi identity is a single scalar equation. One can exploit the 
vector calculus and the differential forms in three dimensions to have a more 
transparent understanding of Hamilton's equations as well as the Jacobi iden- 
tity. Using the isomorphism 

J,=e,,kft^'' = 1,2,3 (2) 

between skew-symmetric matrices and (pseudo)-vectors defined by the com- 
pletely antisymmetric Levi-Civita tensor Sijk, we can write the Hamilton's equa- 
tions and the Jacobi identity as 

V = J xVH J • (V X J) 0, (3) 

respectively. In this form the Jacobi identity is equivalent to the Frobenius 
integr ability condition J AdJ — for the one form J — Jidx^ . It is the condition 
for J to define a foliation of codimension one in three dimensional space |17 ) ,|28 ) - 
[30] . A distinguished property of Poisson structures in three dimensions is the 
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invariance of the Jacobi identity under the multiphcation of Poisson vector J (x) 
by an arbitrary but non-zero factor. The identities 

J ■ V = 0, Vff ■ V = (4) 

follows directly from the Hamilton's equations in ([3]). The second equation in 
(HI is the expression for the conservation of Hamiltonian function. A three 
dimensional vector v (x) is said to be bi-Hamiltonian if there exist two different 
compatible Hamiltonian structures [5].|31). In the notation of equation this 
implies 

V = Ji X VH2 = J2 X VHi (5) 

for the dynamical equations. The compatibility condition for Ji and J2 is 
defined by the Jacobi identity for the Poisson pencil Ji + CJ2 for arbitrary 
constant c. 



3 Frenet-Serret Frame 

Let (t, n, b) denote the Frenet-Serret frame associated with a differentiable curve 
t — i> x{t) in some domain of the three dimensional space R'^. Throughout, 
V — {dx,dy,dz) wiU denote the usual gradient operator in local Cartesian co- 
ordinates. Given a vector field v, the unit tangent vector t, the unit normal n, 
and the unit bi-normal b can be constructed as 

tW-Flii "W-pIl^ b(x)=t(x)xn(x) (6) 

and they form a right-handed orthonormal frame except at points x correspond- 
ing to equilibrium solutions and, those vector fields v satisfying the condition 
imposed by t x (V x t) = 0. This condition excludes essentially the flows with 
constant unit tangent and the points x at which the unit normal n (hence 
the bi-normal b) have zeros. That is, the cases the Frenet-Serret frame is not 
well-defined. To avoid this we may assume that 

(Vxt)7^A(x)t (7) 

for arbitrary nonzero function A (x). That is, we exclude the dynamical systems 
whose unit tangent vectors are the eigenvectors of the curl operator j32j , |33) . 
We introduce the directional derivatives along the triad (t, n, b) as 

d,=t-V 9„ = n • V 5b = b • V (8) 

so that the variables {s,n,b) are the coordinates associated with the Frenet- 
Serret frame based at the point x. Assuming the Cartesian coordinates are 
functions x = x(s,n.6) of the Frenet-Serret coordinates one concludes easily 
that the Jacobian determinant is non-zero and hence the inverse transformation 
s = s(x), n = n{x.), b = fe(x) exists locally, that is, in a sufficiently small 
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neighborhood of a given point xq G M^. These functions may be obtained by 
integrating the one-forms 

T — t ■ dx, r] = n ■ dx, /3 = b • dx (9) 

the last two of which implies n ^constant and b ^constant when restricted to the 
curve x(i). By inverting equations ([5]) we get the expression V — tds+ndn+hdb 
for the Cartesian gradient in the Frenet-Serret frame. We finally note that 
Eq.® is not the only way to construct Frenet-Serret frame, it rather shows 
the existence under the condition in Eq. ^ . For the example of the Darboux- 
Halphen system, we will use its symmetries to construct the moving frame. 

4 Jacobi Identity in Frenet-Serret Frame 

It follows from the identity J • v = that the Poisson vector J has no component 
along the unit tangent vector t. Hence, we set 

J = An + Bb (10) 

for unknown functions A (x) and B (x) satisfying ^ 0. Assuming A^Q 

and defining the function /i — B/A the Jacobi identity for J =A(n-|-/ib) reduces 
to the Riccati equation 

dsn = nn + ^mnh + A*'^b (n) 

in the arclenght variable s. Here, we define the helicity densities 

■Hn = n- (V X n) , Hnb = n-V x b + b-V x n, U\,^ b- (V x b) 

associated with the Frenet-Serret triad [23 • These functions are the coefficients 
of the Cartesian volume dx l\dy K dz in the three-forms 

= 77 A drj, rinh ^r]Adl3 + /3A drj, fib ^ /3 A d(3. (12) 

We refer to [31]- [23 for physical and geometric interpretations of these quanti- 
ties in the context of hydrodynamics. The Riccati equation (fTTj) is equivalent 
to a linear second order equation and hence possesses two linearly indepen- 
dent solutions leading to two Poisson vectors for the dynamical systems under 
consideration. The Hamiltonian form of dynamical equations implies that the 
Poisson vectors obtained from solutions of Riccati equation are always compat- 
ible. Thus, we conclude that 

Proposition 1 All dynamical systems in three dimensions subjected to the con- 
dition ^ posses two compatible Poisson vectors. 
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5 Local Form of Hamiltonian Structures 



By constant solution of Eq. (|lip we mean a function of normal coordinates n 
and h. In particular, such solutions of the Riccati equation include the case 

H„-0, = 0, H„b=0 

the first two of which are the conditions for the normal vectors n and b to satisfy 
the Jacobi identity or Frobenius integrability criterion. The last equation turns 
out to be the compatibility condition for the normal vectors to form a Poisson 
pencil and thereby to define bi-Hamiltonian structure. Integrability criterion 
implies the existence of one-forms ^ and <; such that dr] = S, Arj and d/3 = <; A /3. 
From the compatibility condition we see that ^ = and the common integrating 
factor ^ satisfies ij A = 0, /3 A = 0. Thus, we have d^ = gr] A fi for some 
arbitrary function g and this closes the algebra of one-forms. The integrating 
factor ^ is itself integrable, that is, ^ A = when g — oi equivalently 
= 0. By Poincare lemma, the integrating factor is the differential of a 
function and this makes possible to integrate ry and /? explicitly in the given 
coordinate representation. These results can be drawn from the requirement 
that the one-form 

^global = ?7 + /?£ + fde 

with £ being a parameter, is integrable Tgiotai A dT global = to all orders in 
£ [17] . Eventually, all the known examples of bi-Hamiltonian systems in three 
dimensions are included in this case. The normal coordinates n and b obtained 
from explicit integration of the one-forms rj and (3 represent the global con- 
served quantities. They appear arbitrarily in the Riccati equation and hence in 
the Poisson structures. Conversely, if we are given a bi-Hamiltonian dynamical 
system of the form v =ipVHixVH2 as, for example, in |6]-[26], the representa- 
tion in the Frenet-Serret frame requires t =ip \\ v ||~^ Vi?ixVi?2 which implies 
the orthogonality of the unit tangent vector to the gradients of Hamiltonian 
functions. The normal vectors n and b can then be identified via orthonormal- 
ization procedure applied to the linearly independent vectors Vi?i and \7H2 
taking also the constraint || ViJixViJ2 || = || v || /ip into account. Thus, the 
normal vectors n and b defining the bi-Hamiltonian structure corresponds, in 
the globally integrable case, to the gradients of Hamiltonian functions defining 
Poisson vectors. Local structure arises when we are not able to integrate the 
conserved functions explicitly. Before discussing this case, we remark that, as 
far as the Frenet-Serret frame is constructable, one can still use the normal 
vectors to cast the dynamical equations into formal bi-Hamiltonian form. To 
summarize, we have 

Proposition 2 The manifestly bi-Hamiltonian equation t = n x b is the local 
form of the bi-Hamiltonian structure in three dimensions. Using an appropriate 
volume three-form v, this equation can be written as 

iQV^rjAP (13) 
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where the one-forms 77 and /3 are defined by Eq.(^. 

In the case that the solution fi of the Riccati equation is not a constant, at 
least one of the three-forms in Eg. (1121) is non-zero. Thus, they arise as obstruc- 
tions to extent globally the local bi-Haniiltonian structure. To our knowledge, 
there is only one example in the literature for which we can construct explic- 
itly the obstruction to existence of global Hamiltonian structure. This is the 
Darboux-Halphen system which was obtained from the theory of surfaces more 
than a hundred years ago and resurrected in theoretical physics around eighties. 



6 The Darboux-Halphen System 

In 1878, Darboux 38 obtained, in his investigation of family of second degree 
surfaces in R'^, the system of equations 

d{x + y) d{y + z) d{x + z) 

; = xy , = yz , = xz. (14) 

dt ^ ' dt ^ ' dt ^ ' 

Soon after, Halphen [39], [40] has given the time-dependent transformations 

with {x^,x'^,x^) = {x,y,z), which leave Eq. p^ invariant. He has also ob- 
tained the solution of the system. Starting from 1979, growing attractions on 
the Darboux-Halphen equations come from the observations that various 
important models of theoretical physics admit reductions to equations ([T^ or 
its generalizations. It appeared in the analysis of SO(3) invariant anti-self-dual 
Einstein metrics in gravitation theory j41j , in connection with the moduli space 
of two-monopole problem [42], [43], the complex Bianchi IX cosmological models 
and the reductions of self-dual Yang-Mills fields the WDVV equations of 
topological field theories [JS] and, renormalization group flow of WZW theory 
[46] . It manifests several interesting properties concerning integrability. Namely, 
the Darboux-Halphen system is an algebraically non-integrable system |47).|48) 
whose general solution can be expressed in terms of elliptic integrals |39) . 

6.1 Halphen's symmetries 

We let the Darboux-Halphen system be represented in a local coordinate 
system x = (x,y, z) by the vector field 

d d d 

V = {yz - xy - xz)— + {xz ~ xy - yz)— + {xy - xz - yz)— (16) 
ox ay oz 



for which the corresponding differential equations are equivalent to Eq. (|14p pro- 
vided we perform the replacement t H> —t/2. For notational convenience, we 
also define the autonomous vector fields 
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9 d d . d d d 

" = 2(a;-^ + + ^T^j , — + — + — (17) 

ox Oy oz Ox Oy Oz 

and observe that, they satisfy the Lie bracket relations 

[u, v] = 2v, [u, w] — —2w, [v, w] — u 

of the three-dimensional Lie algebra s[(2). We call the set {v,u,w) of vector 
fields having weights (—1, 0, 1), respectively, an s[(2)— triple. The generators of 
the Halphen's symmetries are the vector fields 

Vt = ^, Ut = -2t^^+u, Wt = -f^^~w + tu (18) 

on the time-extended space / x A/ C M x R"^. These vector fields satisfy the 
conditions 

[|+.,^.]=0, [|+«,t/.] = -2(| + .), = (19) 

to be infinitesimal time-dependent geometric symmetries of the Darboux-Halphen 
vector field v. Eq.([T9]) guarantee that if a function 7J is a time-dependent con- 
served quantity for v, then X{H) with X being any of the vectors in Eq. p^ . 
is also a conserved quantity. Geometric consequences of the time-dependent 
symmetries are discussed in [33]. Below, we will use these symmetries to show 
the local existence of time-independent conserved quantity. This will arise as a 
result of a locally integrable one-form. However, we will not be able to integrate 
it explicitly because the integrating factor turns out to be not integrable in the 
sense of Frobenius. 

Having the sl(2)— triple {v,u,w) intrinsically associated with the Darboux- 
Halphen phase space, one can define a Frenet-Serret frame by orthonormalizing 
these vector fields and proceed as described in previous sections. Namely, one 
can introduce a Poisson vector on the plane of normal vectors, obtain a Riccati 
equation from the Jacobi identity and investigate its solutions. It is obvious from 
the explicit form of vector fields that the coefficients of the three- forms in Eq. ([T^ 
will be non- vanishing. In spite of the explicit expressions, this computational 
way is quite cumbersome and will not be much useful unless one guarantees 
the explicit integration of the one-forms defining the moving coordinates. That 
means, one presumes existence of global bi-Hamiltonian structure. Instead, we 
will first cast the Darboux-Halphen system into the form of a local Hamiltonian 
structure as described by Eq. ([T3l) and see if we can integrate the one-forms 
involved, that is, if we can extend globally the local structure. For this purpose, 
we will proceed by exploiting the Halphen symmetries further. 

A direct computation shows that Eq. (fT9|) are equivalent to 

n o o 

[- + V, Vm] - + Um] = + ^m] = (20) 
where the time-dependent vector fields 

Vm = -V , Um = u + 2tv , Wm = -w + tu + t^v (21) 
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on M C are the unique characteristic (or evolutionary) forms [3] oiVt,Ut, Wt 
along the Darboux-Halphen field v. One can check easily the following result. 

Propositions Each of the sets {v,u,w), {Vt,Ut,Wt) and {Vm,Um,Wm) of 
vector fields form sl{2)~ triples. 

It follows from Eas. (fT8| or ((2T|) that the function p defined by 

p-i = dct((l,v),V,,Ut,Wt) 

= det((l, v), (1, 0), {-2t, u), ( - -w + tu)) 

= det((l,v),Vj^„UM,WM) 

= det((l,v),-v,u + 2tv, - w + tu + t^v) 

= — (u X v) • w = —A{x — y){y — z){z — x) 

is the so called last multiplier [SU]- Here, the boldface letters u, v,w are the 
components of u,v,w, respectively. 

Proposition 4 j^(x) — pdx A dy A dz is the invariant volume element on the 
Darboux-Halphen phase space. 

In other words, with respect to the volume v all three vector fields u,v,w 
are divergence- free. This will be seen easily when we complete the algebra of 
dual one-forms. Let (/3,a,7) be one-forms dual to the s[(2)— triple {v,u,w) so 
that they satisfy 

Zi,/3 = iua = iu,7 = 1 

with the rest of the pairings being zero. Using the components of the rep- 
resentative vector fields in Eqs.([T6)) and (fT7|) of the chosen basis we find the 
expressions 

/3 = pu X w ■ dx, a = pw x v • dx, 7 = pv x u • dx 

from which the orthonormality relations follow easily. For these dual basis one- 
forms, corresponding to the Lie bracket relations of (f , u, w), the s((2)— Maurer- 
Cartan equations 

dl3 = -2a A/3, da = j A l3 , d-f = 2aA-f (22) 

can be obtained using the invariant definition 

2duj{u,v) = uiiyij) — v(iu(-L>) — i[u,v]^ 

of the exterior derivative on a one-form uj and the Lie algebra relations. 

Proposition 5 The Darboux-Halphen dynamical system can be written as 

iyp = a Aj = d'j/2 

which is the form of local Hamiltonian structure. 



8 



Similar relations for the other vector fields of the s[(2)— triple are 

*«A* = 7 A /3 = da, i^l^ = /? A a = dj3/2 

and these relations verify also the invariance of the volume form under the flows 
of s[(2)— triple. It also follows that 

,9 A d/3 = 0, a ^da^ P^^i^, 7 A ^7 = 

so that the one-forms /3 and 7 are integrable in the sense of Frobenius. To this 
end, we recall the integration of one-forms by homotopy formula. Given a one- 
form uj = u}i{x)dx\ let H{x) — x^uJi{sx)ds. Then, dH = a; [3]. However, for 
each of the one-forms /3 and 7 the integrands pu x w ■ x and pv x u • x vanish 
identically and one cannot apply the homotopy formula for their integration. 

6.2 Non-integrable integrating factors 

The algebraic structure on the solution space of the Darboux-Halphen system 
can be obtained, as in the case of globally bi-Hamiltonian structures, from a 
locally integrable one-form involving an arbitrary parameter, this time to second 
order. More precisely, we consider the one form 

Tzocai — 13 + {2am — dm)e — rri'^e^ + (m + ne)de 

and require this to satisfy the local integrability condition Tiocai A dTiocai = 
to all orders in the parameter e. The resulting equations, beyond immediate 
consequences of the Maurer-Cartan equations, include relations for integration 
of locally integrable one-forms. Referring to |T7] for details of this computation, 
we find that the one-forms 

e-^^'^j, e'^'^P (23) 

are closed. Provided these are not elements of the first cohomology class and, by 
Poincare lemma, they must be locally exact, that is, e~^-^"7 = dHi for some 
function Hi and similarly for the other. Since, i„7 = by duality relations, 
this function is a conserved quantity of the Darboux-Halphen system and, as 
we pointed out earlier, the Halphen symmetries, when acted upon Hi produces 
new conserved quantities two of which will be sufficient to establish the bi- 
Hamiltonian structure. For a given parametrized path the integral of the one- 
form a can be computed locally and hence the local conserved quantity along 
the given path can be obtained. This structure remains only local because the 
integrating factor a, for explicit integration of Hi, is non-integrable 

a A da 7^ 

and from the sl(2)— Maurer-Cartan equations, we see that this condition is 
equivalent to the non- vanishing of the Godbillon-Vey three- form |29j.|30j 

a A;9 A7 7^ . 
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